Durum wheat genome highlights past domestication signatures and future improvement targets.

The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

Seed Dormancy Involves a Transcriptional Program That Supports Early Plastid Functionality during Imbibition.

Red rice fully dormant seeds do not germinate even under favorable germination conditions. In several species, including rice, seed dormancy can be removed by dry-afterripening (warm storage); thus, dormant and non-dormant seeds can be compared for the same genotype. A weedy (red) rice genotype with strong dormancy was used for mRNA expression profiling, by RNA-Seq, of dormant and non-dormant dehulled caryopses (here addressed as seeds) at two temperatures (30 °C and 10 °C) and two durations of incubation in water (8 h and 8 days). Aim of the study was to highlight the differences in the transcriptome of dormant and non-dormant imbibed seeds. Transcript data suggested important differences between these seeds (at least, as inferred by expression-based metabolism reconstruction): dry-afterripening seems to impose a respiratory impairment onto non-dormant seeds, thus glycolysis is deduced to be preferentially directed to alcoholic fermentation in non-dormant seeds but to alanine production in dormant ones; phosphoenolpyruvate carboxykinase, pyruvate phosphate dikinase and alanine aminotransferase pathways appear to have an important gluconeogenetic role associated with the restoration of plastid functions in the dormant seed following imbibition; correspondingly, co-expression analysis pointed out a commitment to guarantee plastid functionality in dormant seeds. At 8 h of imbibition, as inferred by gene expression, dormant seeds appear to preferentially use carbon and nitrogen resources for biosynthetic processes in the plastid, including starch and proanthocyanidins accumulation. Chromatin modification appears to be a possible mechanism involved in the transition from dormancy to germination. Non-dormant seeds show higher expression of genes related to cell wall modification, suggesting they prepare for acrospire/radicle elongation.

Comparative Transcriptome Profiles of Near-Isogenic Hexaploid Wheat Lines Differing for Effective Alleles at the 2DL FHB Resistance QTL.

Fusarium head blight (FHB), caused by the fungus Fusarium graminearum, represents one of the major wheat diseases worldwide, determining severe yield losses and reduction of grain quality due to the accumulation of mycotoxins. The molecular response associated with the wheat 2DL FHB resistance QTL was mined through a comprehensive transcriptomic analysis of the early response to F. graminearum infection, at 3 days post-inoculation, in spikelets and rachis. The analyses were conducted on two near isogenic lines (NILs) differing for the presence of the 2DL QTL (2-2618, resistant 2DL+ and 2-2890, susceptible null). The general response to fungal infection in terms of mRNAs accumulation trend was similar in both NILs, even though involving an higher number of DEGs in the susceptible NIL, and included down-regulation of the primary and energy metabolism, up-regulation of enzymes implicated in lignin and phenylpropanoid biosynthesis, activation of hormons biosynthesis and signal transduction pathways and genes involved in redox homeostasis and transcriptional regulation. The search for candidate genes with expression profiles associated with the 2DL QTL for FHB resistance led to the discovery of processes differentially modulated in the R and S NILs related to cell wall metabolism, sugar and JA signaling, signal reception and transduction, regulation of the redox status and transcription factors. Wheat FHB response-related miRNAs differentially regulated were also identified as putatively implicated in the superoxide dismutase activities and affecting genes regulating responses to biotic/abiotic stresses and auxin signaling. Altered gene expression was also observed for fungal non-codingRNAs. The putative targets of two of these were represented by the wheat gene WIR1A, involved in resistance response, and a gene encoding a jacalin-related lectin protein, which participate in biotic and abiotic stress response, supporting the presence of a cross-talk between the plant and the fungus.

Wild emmer genome architecture and diversity elucidate wheat evolution and domestication.

Wheat (Triticum spp.) is one of the founder crops that likely drove the Neolithic transition to sedentary agrarian societies in the Fertile Crescent more than 10,000 years ago. Identifying genetic modifications underlying wheat's domestication requires knowledge about the genome of its allo-tetraploid progenitor, wild emmer (T. turgidum ssp. dicoccoides). We report a 10.1-gigabase assembly of the 14 chromosomes of wild tetraploid wheat, as well as analyses of gene content, genome architecture, and genetic diversity. With this fully assembled polyploid wheat genome, we identified the causal mutations in Brittle Rachis 1 (TtBtr1) genes controlling shattering, a key domestication trait. A study of genomic diversity among wild and domesticated accessions revealed genomic regions bearing the signature of selection under domestication. This reference assembly will serve as a resource for accelerating the genome-assisted improvement of modern wheat varieties.

microRNAs differentially modulated in response to heat and drought stress in durum wheat cultivars with contrasting water use efficiency.

Plant stress response is a complex molecular process based on transcriptional and posttranscriptional regulation of many stress-related genes. microRNAs are the best-studied class of small RNAs known to play key regulatory roles in plant response to stress, besides being involved in plant development and organogenesis. We analyzed the leaf miRNAome of two durum wheat cultivars (Cappelli and Ofanto) characterized by a contrasting water use efficiency, exposed to heat stress, and mild and severe drought stress. On the whole, we identified 98 miRNA highly similar to previously known miRNAs and grouped in 47 MIR families, as well as 85 novel candidate miRNA, putatively wheat specific. A total of 80 known and novel miRNA precursors were found differentially expressed between the two cultivars or modulated by stress and many of them showed a cultivar-specific expression profile. Interestingly, most in silico predicted targets of the miRNAs coming from the differentially expressed precursors have been experimentally linked in other species to mechanisms controlling stomatal movement, a finding in agreement with previous results showing that Cappelli has a lower stomatal conductance than Ofanto. Selected miRNAs were validated through a standardized and reliable stem-loop qRT-PCR procedure.

The genome sequence of the outbreeding globe artichoke constructed de novo incorporating a phase-aware low-pass sequencing strategy of F1 progeny.

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F1 individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.

The Citrus clementina putative allergens: from proteomic analysis to structural features.

Several allergens have been identified and characterized in the genus Citrus, which belongs to the germin-like proteins (GPLs), profilins, and non-specific lipid transfer proteins (nsLTPs). In this work, in silico sequence analysis, protein purification, mass spectrometry identification, and the spectral counting method were integrated to identify new putative allergens of Citrus clementina and their expression level in the fruit peel. The in silico analysis revealed fifteen new sequences belonging to GLPs (Cit cl 1), and two more belonging to nsLTPs (Cit cl 3). No other new sequences were found as regards profilins (Cit cl 2). Each putative allergen from fruit peel was obtained using different protein extraction methods, and the protein sequences of the putative allergens were identified by means of LTQ-Orbitrap XL mass spectrometer. The spectral counting strategy revealed that Cit cl 1 had a higher expression level than Cit cl 2 and Cit cl 3. To predict the quaternary structure and deduced function of Cit cl 1, its primary sequence was used as a template to search a homologous protein structure in the RCSB PDB Database, getting high correspondence with the oxalate oxidase protein in barley.

A Survey of MicroRNA Length Variants Contributing to miRNome Complexity in Peach (Prunus Persica L.).

MicroRNAs (miRNAs) are short non-coding RNA molecules produced from hairpin structures and involved in gene expression regulation with major roles in plant development and stress response. Although each annotated miRNA in miRBase (www.mirbase.org) is a single defined sequence with no further details on possible variable sequence length, isomiRs - namely the population of variants of miRNAs coming from the same precursors - have been identified in several species and could represent a way of broadening the regulatory network of the cell. Next-gen-based sequencing makes it possible to comprehensively and accurately assess the entire miRNA repertoire including isomiRs. The aim of this work was to survey the complexity of the peach miRNome by carrying out Illumina high-throughput sequencing of miRNAs in three replicates of five biological samples arising from a set of different peach organs and/or phenological stages. Three hundred-ninety-two isomiRs (miRNA and miRNA*-related) corresponding to 26 putative miRNA coding loci, have been highlighted by mirDeep-P and analyzed. The presence of the same isomiRs in different biological replicates of a sample and in different tissues demonstrates that the generation of most of the detected isomiRs is not random. The degree of mature sequence heterogeneity is very different for each individual locus. Results obtained in the present work can thus contribute to a deeper view of the miRNome complexity and to better explore the mechanism of action of these tiny regulators.

On the complexity of miRNA-mediated regulation in plants: novel insights into the genomic organization of plant miRNAs.

MicroRNAs (miRNAs) are endogenous small non-coding RNAs of about 20-24 nt, known to play key roles in post-transcriptional gene regulation, that can be coded either by intergenic or intragenic loci. Intragenic (exonic and intronic) miRNAs can exert a role in the transcriptional regulation and RNA processing of their host gene. Moreover, the possibility that the biogenesis of exonic miRNAs could destabilize the corresponding protein-coding transcript and reduce protein synthesis makes their characterization very intriguing and suggests a possible novel mechanism of post-transcriptional regulation of gene expression. This work was designed to carry out the computational identification of putative exonic miRNAs in 30 plant species and the analysis of possible mechanisms involved in their regulation. The results obtained represent a useful starting point for future studies on the complex networks involved in microRNA-mediated gene regulation in plants.

First survey of the wheat chromosome 5A composition through a next generation sequencing approach.

Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.

A computational-based update on microRNAs and their targets in barley (Hordeum vulgare L.).

BACKGROUND: Many plant species have been investigated in the last years for the identification and characterization of the corresponding miRNAs, nevertheless extensive studies are not yet available on barley (at the time of this writing). To extend and to update information on miRNAs and their targets in barley and to identify candidate polymorphisms at miRNA target sites, the features of previously known plant miRNAs have been used to systematically search for barley miRNA homologues and targets in the publicly available ESTs database. Matching sequences have then been related to Unigene clusters on which most of this study was based. RESULTS: One hundred-fifty-six microRNA mature sequences belonging to 50 miRNA families have been found to significantly match at least one EST sequence in barley. As expected on the basis of phylogenetic relations, miRNAs putatively orthologous to those of Triticum are significantly over-represented inside the set of identified barley microRNA mature sequences. Many previously known and several putatively new miRNA/target pairs have been identified. When the predicted microRNA targets were grouped into functional categories, biological processes previously known to be regulated by miRNAs, such as development and response to biotic and abiotic stress, have been highlighted and most of the target molecular functions were related to transcription regulation. Candidate microRNA coding genes have been reported and genetic variation (SNPs/indels) both in functional regions of putative miRNAs (mature sequence) and at miRNA target sites has been found. CONCLUSIONS: This study has provided an update of the information on barley miRNAs and their targets representing a foundation for future studies. Many of previously known plant microRNAs have homologues in barley with expected important roles during development, nutrient deprivation, biotic and abiotic stress response and other important physiological processes. Putative polymorphisms at miRNA target sites have been identified and they can represent an interesting source for the identification of functional genetic variability.

Genetic variants of HvCbf14 are statistically associated with frost tolerance in a European germplasm collection of Hordeum vulgare.

Two quantitative trait loci (Fr-H1 and Fr-H2) for frost tolerance (FT) have been discovered on the long arm of chromosome 5H in barley. Two tightly linked groups of CBF genes, known to play a key role in the FT regulatory network in A. thaliana, have been found to co- segregate with Fr-H2. Here, we investigate the allelic variations of four barley CBF genes (HvCbf3, HvCbf6, HvCbf9 and HvCbf14) in a panel of European cultivars, landraces and H. spontaneum accessions. In the cultivars a reduction of nucleotide and haplotype diversities in CBFs compared with the landraces and the wild ancestor H. spontaneum, was evident. In particular, in cultivars the loss of HvCbf9 genetic variants was higher compared to other sequences. In order to verify if the pattern of CBF genetic variants correlated with the level of FT, an association procedure was adopted. The pairwise analysis of linkage disequilibrium (LD) among the genetic variants in four CBF genes was computed to evaluate the resolution of the association procedure. The pairwise plotting revealed a low level of LD in cultivated varieties, despite the tight physical linkage of CBF genes analysed. A structured association procedure based on a general liner model was implemented, including the variants in CBFs,of Vrn-H1, and of two reference genes not involved in FT (alpha-Amy1 and Gapdh) and considering the phenotypic data for FT. Association analysis recovered two nucleotide variants of HvCbf14 and one nucleotide variant of Vrn-H1 as statistically associated to FT.

A combined strategy of "in silico" transcriptome analysis and web search engine optimization allows an agile identification of reference genes suitable for normalization in gene expression studies.

Traditionally housekeeping genes have been employed as endogenous reference (internal control) genes for normalization in gene expression studies. Since the utilization of single housekeepers cannot assure an unbiased result, new normalization methods involving multiple housekeeping genes and normalizing using their mean expression have been recently proposed. Moreover, since a gold standard gene suitable for every experimental condition does not exist, it is also necessary to validate the expression stability of every putative control gene on the specific requirements of the planned experiment. As a consequence, finding a good set of reference genes is for sure a non-trivial problem requiring quite a lot of lab-based experimental testing. In this work we identified novel candidate barley reference genes suitable for normalization in gene expression studies. An advanced web search approach aimed to collect, from publicly available web resources, the most interesting information regarding the expression profiling of candidate housekeepers on a specific experimental basis has been set up and applied, as an example, on stress conditions. A complementary lab-based analysis has been carried out to verify the expression profile of the selected genes in different tissues and during heat shock response. This combined dry/wet approach can be applied to any species and physiological condition of interest and can be considered very helpful to identify putative reference genes to be shortlisted every time a new experimental design has to be set up.

Chromosome regions and stress-related sequences involved in resistance to abiotic stress in Triticeae.

Drought, low temperature and salinity are the most important abiotic stress factors limiting crop productivity. A genomic map of major loci and QTLs affecting stress tolerance in Triticeae identified the crucial role of the group 5 chromosomes, where the highest concentration of QTLs and major loci controlling plant's adaptation to the environment (heading date, frost and salt tolerance) has been found. In addition, a conserved region with a major role in drought tolerance has been localized to the group 7 chromosomes. Extensive molecular biological studies have led to the cloning of many stress-related genes and responsive elements. The expression of some stress-related genes was shown to be linked to stress-tolerant QTLs, suggesting that these genes may represent the molecular basis of stress tolerance. The development of suitable genetic tools will allow the role of stress-related sequences and their relationship with stress-tolerant loci to be established in the near future.